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We investigate the stability of the Bose-Einstein condensate (BEC) for 
the case of atoms with negative scattering lengths at zero temperature using 
the Ginzburg-Pitaevskii-Gross (GPG) stationary theory. We have found a 
new exact equation for determining the upper bound of the critical numbers 
N cr of atoms for a metastable state to exist. Our calculated value of N cr for 
Bose-Einstein condensation of lithium atoms based on our new equation is 
in agreement with those observed in a recent experiment. 
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1. Introduction 

The concept of the Bose-Einstein condensation (BEC) [1] has been known 
for 73 years, and has been used to describe all physical scales, including liq- 
uid 4 He, exitons in semiconductors, pions and kaons in dense nuclear matter 
(neutron stars, supernovae), and elementary particles [2]. It is only a few 
years ago that the BEC phenomenon was observed directly in dilute vapors 
of alkali atoms, such as rubidium [3], lithium [4,5], and sodium [6], confined 
in magnetic trap and cooled down to nanokelvin temperatures. At such low 
temperatures, the thermal de Brogle wavelength increase to submicrometer 
dimensions. When the scattering length is negative, then nonlinear interac- 
tion between atoms is attractive and it has been claimed that the BEC in 
free space is impossible [7], because the attraction makes the system tend 
to an ever denser phase. For 7 Li, the s-wave scattering length a = (-14.5 ± 
0.4) A [8]. For bosons trapped in an external potential, there may exist a 
metastable BEC state under certain conditions for a number of atoms below 
the critical value N cr [9-11]. Stability of the BEC in a trap as a function 
of nonlinearity has been studied numerically for T = [12 - 14]. It was 
demonstrated [15, 16] that global solutions of the Ginzburg-Pitaevskii-Gross 
equation (GPG) [17] for trapped atoms with negative energies diverge in a 
finite time. When the system collapses, a more accurate theory is required 
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in order to include short range effects. 

In this paper, we report our theoretical studies of the stability of the 
BEC for the case with negative scattering length at zero temperature using 
the GPG stationary theory [17]. We have found, for the first time, an exact 
equation for determination of the upper bound of the N cr for a metastable 
state to exist. 

2. A necessary condition for a local minimum 

In the mean-field approximation, the ground state energy of the system 
is given by the GPG energy functional [17] 

Jty)=T + V + ^V H , (1) 



with 



J(V>) = E/N, (2) 

h 2 



J d 3 r^*Vilj, (4) 



T 

V 

v B =l m 4 d\ (5) 



and 



4irh 2 a 

9o = , 6 

m 

where a is the s-wave scattering length, N is the number of particles in the 
BEC, V is an external trapping potential, and ip is the condensate ground 
state wave function normalized as 



/ 



|Vf d 3 r = 1. (7) 
For a harmonic oscillator trap, we have 

V ^ = Vl ( X ) + V 2 (y) +V ^ Z ) ) ( 8 ) 

where 
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Using 

ViiKn) = xfyin), (9) 

we can obtain an extremum condition (see Appendix) 

a N - 

2% - 1% + y -^-V H = 0, (10) 

f^f^m^ (id 
Vi = ^f dhri^i^r^m, (12) 

and a necessary local minimum condition 

A > 0, (13) 
where A is matrix with matrix elements 

Aij = (2Ti + m)5 %3 + (1 - Sij)^-V H . (14) 

(See Appendix for details.) 

Eqs. (10-12) represent well-known virial theorem equations [12], but the 
conditions (13) was not discussed in literature, to the best of our knowledge. 
We note, that if the conditions (13) are not fulfilled, local minimum disap- 
pears and solution will be not stable. In the case a < 0, we can rewrite (14) 

as 

Aij ^ ^ VH [g i 5 ij + 25 ij -l], (15) 

where 

9i = | | K jl-r ( 16 ) 

\9o\NV H 

Matrix (28y — 1) is not positive (its eigenvaliues are -1, 2, 2) hence in case 
Vi — V 2 — V 3 — (free space) matrix A, is not positive and there is no stable 
solution. 

3. Upper bounds for N cr . 
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We consider first an isotropic trap case with 



~ ~ - 16 V 

9x = 9y = 9z = ^\^NV-H (17) 



where 



V = % + V y + V z = 3V X . (18) 

For this case the local minimum conditions given by Eq. (13) can be written 

as 

If > 1. (19) 

3\g \NV H 1 ; 

In terms of N, which is solution of the following equation 

3 \g \V n 

we can state that if N > N, the BEC is not stable. Therefore N is an upper 
bound for N cr , 

N cr < N. (21) 

We note, that Eq.(21) is also valid for anisotropic cases, but the bound given 
by Eq. (20) in this case can be improved. To show this let us consider an 
anisotropic harmonics trap 

V(r) = V ± (f) + V z (f), (22) 

here 

^(0 = |4 o) V + y 2 ), (23) 

V,tf) = y^V. (24) 
Using Eqs. (22-24) we can rewrite the local minimum conditions (13) as 

rW" + (^nH 2 ^ K - 1 > 0, (25) 
\g \NV H \go\NV H 



where 



V ± = ^i o)2 / r^(x 2 + y 2 )^d s f, (26) 
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V z = ^J^ 2 J ^*(f)z 2 ^(f)d 3 f. (27) 
Setting the left-hand side of Eq. (25) to zero we obtain the following equation 

iV = ^[l + (l + ^) 1/2 ] (28) 
\9o\Vh V± 

for upper bound N 

N cr < N. (29) 

It should be pointed out that our formulas (28) and (29) are valid for a 
general anisotropic trap with 

VI = ^ J V>*(f)(c4 o) V + 4 o)2 y 2 )i>(r)dr. (30) 

To the best of our knowledge, these formulas have never been discussed 
previously in the literature. 

4. Application to lithium atoms: numerical example 

We assume the trial wave function for ip in Eqs. (10-28) to be [14, 18] 

iHlf) = w y2 a; l/2 ( Z|)3/4 e -m(.i^ + ^) + ^^)/2,_ (3l) 

Tin 

Using Eq. (31), we rewrite Eqs. (10) and (28) as 

l-25i = 5i(l + 8^A 2 ) 1 / 2 , (32) 

1 - X 2 5 2 z = 5 Z 8±[1 + (1 + 8-^A 2 ) 1/2 ], (33) 

and 

n = ^ ± [l + (l + 8^A 2 ) 1 / 2 ], (34) 
where 8 Z = Wj^/u z , 5± = uj^ /uj±, A = uj^/uj^, and 

n = 2( ^° ))1/2( 2^ )1/2iV|a| - (35) 
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We find that numerical solution of Eqs. (32-34) for < A < 1 can be 
interpolate as 

n = e' ia+ ^ (36) 
with a = 0.490419, f3 = 0.149175. Using Eqs. (36) and (35) we have 

N =(—) 1/2 -2UT- ( 37 ) 

Taking the experimental parameters [5], u^/2tv = 152 Hz, and u>^ /2tt = 
132 Hz, we obtain A = 0.86842 and N cr < N= 1456. This value of N cr is 
consistent with theoretical predictions [12-14] and is in agreement with those 
observed in a recent experiment [5]. 

5. Summary and Conclusions 

We have investigated the stability of the BEC for the case of atoms with 
negative scattering length at zero temperature in a magnetically trapped 
atomic gas. Using a rigorous derivation, we have found a new exact equation 
for determining the upper bound of the critical number of atoms N cr for a 
metastable state to exist. Our calculated value of N cr for BEC of lithium 
is about 1456 atoms, which is consistent with recent experimental measure- 
ments [5] and theoretical predictions [12-14]. 
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Appendix 



In this appendix, we give proving Eqs. (10-14). Let us introduce a trial 
function 

<Mf) = [(1 + Cl )(l + e 2 )(l + e 3 )] 1/2 ^((1 + e 1 )x, (1 + e 2 )y, (1 + e 3 )z), (A.1) 

where tp(x,y,z) is exact ground state solution of the GPG equation. The 
evaluation of J(V't), Eq. (1), can be carried out as 



+ 



where 



J{ip t ) = J(ei,e 2 ,e 3 ) = 



(1 + ^) 2 



7<. 



1 + e 1 ) 2 f x + (1 + e 2 ) 2 f, + (1 + e 3 ) 2 f z + 

1 t/ 2 + (l + e 1 )(l + e 2 )(l + e 3 )^#t^, 



/ 



l + e 3 ) 2 

h 2 



2 

(A.2) 



cfr 



2m 



and 



Vi 



m 



iv^r, 



y /rf 3 r^*(f)(4° )2 r 2 )^iH. 



Since ^(r) is exact ground state solutions, J(ei,e 2 ,e 3 ) has an extremum at 
£i — e 2 = ^3 = 0, which requires that 



<9J(ei,e 2 ,e 3 ) 



Substitution (A.2) into (A.3) yields Eq. (10), 

a N - 
2Ti -2Vi + ^— V H 



= 0. 



0. 



(A3) 



(A4) 



Eq. (A. 4) represents well-known virial theorem equations [12]. 
Now, let us introduce matrix A with matrix elements 



<9 2 J(ei,e 2 ,e 3 ) 



deidej 



- a N - 

(2T i + 6V i )5 ij + (l-5 ij )^V H . 



ei,e 2 ,e 3 =0 

In terms of A we can write a local minimum conditions as 

A > 



(A.5) 



(A6) 
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